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Abstract 

We apply the homotopy group theory in classifying the topological defects in 
atomic spin-1 and spin-2 Bose-Einstein condensates. The nature of the defects de- 
pends crucially on the spin-spin interaction between the atoms. We find the topolog- 
ically stable defects both for spin-1 ferromagnetic and anti-ferromagnetic states, and 
for spin-2 ferromagnetic and cyclic states. With this rigorous approach we clarify the 
previously controversial identification of symmetry groups and order parameter spaces 
for the spin-1 anti-ferromagnetic state, and show that the spin-2 cyclic case provides a 
rare example of a physical system with non-Abelian line defects, like those observed 
in biaxial nematics. We also show the possibility to produce vortices with fractional 
winding numbers of 1/2, 1/3 and their multiples in spinor condensates. 



1 Introduction 

The all-optical trapping of Bose-Einstein condensates (BEC) U |2 ] has opened up a new 
direction in the study of dilute atomic gases, i.e., the spinor condensates with degenerate 
internal degrees of freedom of the hyperfine spin F. For alkali atoms with F = 1, both 
experiments and theories have shown two possible kinds of spin correlations in the atom 
species, namely ferromagnetic (e.g. 87 Rb 00113) or antiferromagnetic (e.g. 23 Na 
With the experimental success of condensing alkali bosons with F > 1 such as 85 Rb Q 
and 133 Cs (SJ, and the unusual stability of the F = 2 state (against spin-exchange) in 
87 Rb [9 1, one expects that defects with much richer structure can be created in the future. 
A remarkable feature here is that both the gauge symmetry U(l) and the spin symmetry 
SO(3) are involved, a situation similar to superfluid 3 He where three different continuous 
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symmetries (orbital, spin and gauge) are broken either independently or in a connected 
fashion fTOllTTl . 

Topological defects and excitations in the spinor BECs have been studied theoretically 
by several groups El |H [121 [131 QH E21 EU [TTI HU EU |20|. Stoof and Khawaja E) 
showed that ferromagnetic condensates have long-lived Skyrmion excitations, which are 
nonsingular but topologically nontrivial pointlike spin textures. Moreover, they also found 
that spin-1 Bose-Einstein antiferromagnets have singular pointlike topological spin tex- 
tures [13 1, which are analogous to the 't Hooft-Polyakov magnetic monopoles in particle 
physics. Coreless vortices were demonstrated to be fhermodynamically stable in ferromag- 
netic F = 1 spinor condensates under rotation ff4l l~j~6l and were phase imprinted in a 
F = 1 sodium condensate experimentally 1211 . Yip flD has performed a systematic study 
on vortex structures and presented several axisymmetric and non-axisymmetric vortices for 
F = 1 antiferromagnetic BEC. Martikainen et al. IIT31 proposed and demonstrated nu- 
merically a method to create monopoles in three dimensional two-component condensates. 
Linear defects were studied by Leonhardt and Volovik lTT8l . who pointed out the existence 
of Alice strings in the condensate of 23 Na. 

Most of the work on this subject is based on the original identification of the order pa- 
rameter spaces by Ho [3]. After the original studies it was also claimed by Zhou that a 
discrete symmetry of Z2 type was missed in the case of antiferromagnetic spin-1 conden- 
sate lfT9l |20 ] and therefore the topological defects would manifest totally different struc- 
tures. In this article we present a rigorous topological study that both solves this spin-1 
controversy, and reveals interesting aspects of spin-2 systems. The phases of spin-2 spinor 
condensates are characterized by a pair of parameters |(F)| and |0| describing the ferro- 
magnetic order and the formation of singlet pairs, respectively l22l l23l l24l l25l . It turns 
out that for the so called cyclic phase the fundamental group that determines the nature of 
possible stable topological defects is non-Abelian. The only known physical example of 
such a system so far has been the biaxial nematic liquid crystal. 

The organization of this paper is as follows: In the following section, we shall review the 
basic physics of the spinor condensate and discuss the possible ground states for hyperfine 
spin F = 1 and F = 2. In Section 3, we give a brief introduction of the homotopy theory of 
the defect classification, taking the nematic liquid crystal and superfluid 3 He as examples. 
We present our calculation of the homotopy groups for spinor condensates in Sections 4-7. 
The non-Abelian fundamental group for the cyclic phase and its indications are discussed 
in detail and the order parameter spaces are easily identified in a correct way following our 
procedure of symmetry breaking. We summarize our results in Section 8. 

2 Spinor Condensate 

Neutral atomic gases can be confined in conventional magnetic traps with the availability of 
hyperfine states being restricted by the requirement that the trapped atoms remain in weak- 
field seeking states. Alkali atoms with a nuclear spin of / = 3/2, such as 87 Rb and 23 Na, 
have three weak-field seeking states at small field. A far-off-resonant optical trap, however, 
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confines atoms regardless of their hyperfine state. Thus, the atomic spin is liberated from 
the requirements of magnetic trapping and becomes a new degree of freedom. In particular, 
all atoms in the lower hyperfine manifold, for example the F = 1 hyperfine manifold 
of sodium, can be stably trapped simultaneously. Such multi-component optically uapped 
condensates are represented by an order parameter which is a vector in hyperfine spin space, 
and are thus called spinor Bose-Einstein condensates. The spin relaxation collisions in 
spinor condensates allow for population exchange among hyperfine states without trap loss. 
Theoretical studies started with the determination of the ground state structure in mean field 
theory for both spin-1 and spin-2 E2 |23] M E3 HD E2 cases. Law et al. I28ll29l 

investigated the spin correlation beyond mean-field limit and the spin-mixing dynamics 
due to the nonlinear interaction in the spinor condensate. The dynamics is sensitive to the 
relative phase and particle number distribution among the individual components of the 
condensate. Ho and Yip l30l later found that the ground state of a spin-1 Bose gas with 
an antiferromagnetic interaction was a fragmented condensate in uniform magnetic fields. 
Zhou [ 19 20 1 showed that the low energy spin dynamics in the system can be mapped into 
an o(n) nonlinear sigma model. The formation of ground state spin domains, metastable 
states and quantum tunneling were observed in a series experiments at MIT 161 I3T1I321I331 . 
The discussions in this paper, however, mainly concern the possible ground states in mean- 
field theory. 

2.1 Spin-1 case 

The ground states of the spinor condensate are determined through the minimization of the 
energy functional with the constraint of the conservation of the atom number and magneti- 
zation [33 1 . An F = 1 spinor Bose-Einstein condensate is described by a three-component 
order parameter ?/>(r) = ^(r)^) = (tp+i, ipo, i J -i) T - I n second quantized notation, the 
Hamiltonian describing a weakly-interacting Bose gas can be obtained from the Gross- 
Pitaevskii theory 

H = |d 3 r|*J(r)(-^ + C/(r))^(r)^ 

+- so *t(r)*t(r)* i (r)* i (r) 

+lg 2 ^(r)^(r)(F a ) ik (F a ) jl i> k (r)^ l (r)\ (1) 

where \&j(r) is the field annihilation operator for an atom with mass m in hyperfine state 
at position r with i = +1,0,-1 and U(r) is the trapping potential. Here the re- 
peated indices are summed. The scattering lengths ao and a 2 characterize collisions be- 
tween atoms through the total spin and 2 channels, respectively, go = a +2a 2 - ls 
interaction strength through the "density" channel, and g 2 = a2 ~ a ° is that through the 
"spin" channel. 
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It is convenient to express the order parameter as ^(r) = -y/n(r)£(r) where n(r) 
is the atomic density and £(r) is a three-component spinor £(r) = (C+i, Co> C-i) T = 
(x + e ie+ , x^e 1 ® , x_e l6l ~) of normalization |£| 2 = 1. Here x and 9 are the amplitudes 
and phases of the components. The spinor determines the average local spin by means of 
(F) = £t(r)F£(r), and F are the usual spin-1 matrices with 



F x 



which obey the commutation relations [F a , F^] = ie a b c F c . We thus obtain the energy func- 
tional 
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J d 3 r (-^-) i> + (U(r) - „) n + \ (g + g 2 (F) 2 ) } 

= J d 3 r (Kq + n 2 g 2 (F) 2 /2) (2) 

where i\~o is the density-dependent part and the chemical potential fi determines the number 
of atoms in the condensate. It is obvious that all spinors related to each other by gauge trans- 
formation e %6 and spin rotations U = e ~ lFza e~ lFyl3 e~ lFzl are energetically degenerate in 
zero external magnetic field, where (a, (3, 7) are the Euler angles. The ground-state spinor is 
determined by minimizing the spin-dependent mean-field interaction energy, n 2 g2 (F) 2 /2. 
There are two distinct states depending on the sign of the interaction parameter g 2 : 

• 92 > (i.e. 02 > ao, e.g. 23 Na): anti-ferromagnetic or polar state as the condensate 
lowers its energy by minimizing its average spin, i.e. by making (F) = 0. The 
ground state spinor is then one of a degenerate set of spinors, the "polar" states, 
corresponding to all possible rotations of the hyperfine state nip = 0, i.e. 



C(r) = e M U 



1 ^— ia 



( 11 / \ 




e %a sin (3 
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• 92 < (i.e. 02 < ao, e.g. 87 Rb): ferromagnetic as the condensate lowers its energy 
by maximizing its average spin, i.e. by making (F) = 1. In this case the ground state 
spinors correspond to all rotations of the hyperfine state rap = 1, i.e. 




e- ia cos 2 I 



C(r) = e*W| |= e ^-T)| V2cosfsinf ] (4) 

e ia sin 2 4 
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2.2 Magnetic Field 



One can tailor the ground state structure with an external magnetic field and the effects of 
field inhomogeneities and quadratic Zeeman shifts modify the spin-dependent interaction 
energy into [ 6 ] 

K spin = (c (F) 2 - p (F z ) + q <F 2 )) n (5) 

where c = ginj2. The linear Zeeman shift p = q^bBz+pq, where g is the Lande g-f actor 
and [is is the Bohr magneton, comes from the field gradient B along the long axis z of the 
condensate, while the last term gives the quadratic Zeeman shift from homogeneous field 
which is always positive for spin-1 condensate in a weak field. Assuming conservation 
of total spin, we have included a Lagrange multiplier po into p. For a system with zero 
total spin, po cancels the linear Zeeman shift due to a homogeneous bias Bo, yielding 
p = 0. Positive (negative) values of p are achieved for condensates with a positive (negative) 
overall spin. The parameters p and q can be related to the individual level shifts by (energies 
in units of the hyperfine splitting Ehfs) 

2p = E + 

2q = E_ + E + - 2E (6) 

where the Zeeman energies E + , E$ and EL of the mp = +1,0,-1 can be expressed 
according to the Breit-Rabi formula l34l as 

1 1 / o 

E+ = ----Vl + x + x 2 

o I 

E = -I-IyiT^ 2 

E- = -^-^x/l-x + x 2 (7) 

with x = gfi B B/E H FS- 

Including the non-diagonal terms of the mean field interaction, we may minimize the 
energy functional 

Kspin/n = c (x 2 + — x 2 ) 2 + 2cxq (x 2 + + x 2 + 2x + x^ COS (f) 

-p(x 2 + -x 2 _) +q{x\ +x 2 _) (8) 

by means of the Lagrange multiplier method subjected to the constraint of normalization 

g = x\ + x 2 + x 2 - 1 = (9) 

where <p = 9 + + 9- — 29q. The solutions to the first derivatives of the Lagrange multiplier 
function X = K sp i n /n — Xg can be classified into the following table of spinors with their 
corresponding energies 
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spmors 



energies 



1 

2 
3 



e 4 ^ 1 (1,0,0) 
e^- 1 (0,0,1) 
e id ° (0,1,0) 



c — p + q 
c+p + q 
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(2c+p+g) 2 



(2c-p+g) 2 



(g 2 +4cg-p 2 ) 
ldcq 2 
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We notice that spinors 4-7 are only well-defined in some specific regions in the p-q 
plane, i.e., the quantities under the square root must be non-negative. For example, spinor 7 
may only exist for q 2 + Acq — p 2 < and q 2 > p 2 , and in addition we must have cj> = or 
7T. The ground state spinors obtained by minimizing the energy functional can be indicated 
in the so-called spin-domain phase diagrams (Figure 1 in ref. 161). For c = 0, the Zeeman 
energy causes the cloud to separate into three pure domains with mp = +1,0,— 1 and with 
boundaries at \p\ = q. For c > 0, a spin domain with mixed mp = ±1 components, i.e., 
spinor 4, appears in the anti-ferromagnetic phase diagram. For c < 0, all three components 
are generally miscible and have no sharp boundaries, which corresponds to spinor 7. 

2.3 Conservation of Magnetization 

Although conservation of the magnetization was included in the above section, it was not 
separately discussed. Consequently the results do not easily apply to systems with fixed 
values of the magnetization M. . The ground state structures as given in [6 1 correspond to the 
actual ground state as realized through an M. non-conserving evaporation process (e.g. in 
the presence of a non-zero l?-field) that serves as a reservoir for condensate magnetization. 
On the other hand, the phase diagram for fixed values of M. was also explicitly discussed 
ll35l . which could physically correspond to experimental ground states (with/without a B- 
field) due to an J\A conserving evaporation process. This requires the introduction of two 
Lagrange multipliers during the minimization subjected to conservation constraints for both 
the atomic number N and magnetization M., which in the mean-field approximation are 
given by 



M 



N = 




(10) 
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We restrict the discussion here to the situation that equation © and 
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(11) 



are satisfied where m = M./N. With the definition x 



+ x_, we can assort the 



possible spinors minimizing the Lagrange multiplier function X = K spin /n — Xg — 5h into 
the following classes (where the energy zero point has been moved to pm): 



spinors 



energies 



1+m 



m, e 



, ,0,e« 



l—m 



i,0) 

(0, e ie ° y/l + m, e 19 - 1 yf^rn) 



cm? + 



i0 + l I Xm+m i6 



J9. 



—cm? + (2c + q)m 
—cm 2 — (2c + q) m 

cm? + g±(x m ) + qx r 



We still have the spinors 1-3 which are the same as in above section, however, they only 
exist for special values m = +1,-1,0, respectively. While spinor 5(6) is confined to the 
positive (negative) values of to, 4 and 7 may exist for the whole region — 1 < m < 1. In 
spinor 7 with three nonzero components, the phase convention remains cj) = or n and the 
minimum is reached when x = x m where x m is determined by 

g' ± {x m ) + q = (12) 

with g±(x) = 2c(l — x) (x ± \J x 2 — m 2 ^j for ferromagnetic(+) or anti-ferromagnetic(— ) 
interaction, respectively. The ground state spinor phase diagram for a homogeneous con- 
densate may be determined in the m-q plane, as indicated for positive to case in Figure 4 of 
ref. l35l . For c = 0, spinor 5 will always dominate except that on the boundary q = we 
have spinor 7. For c < 0, spinor 7 will dominate. For c > 0, a curve q = 2c(l — \/l — m 2 ) 
divides spinors 4 and 7. 



2.4 Spin-2 case 

For 23 Na and 87 Rb with regular hyperfme multiplets, the lower hyperfine state F = 1 has 
lower energy than the upper state F = 2. Experimentally only atoms in the lower hyperfine 
states can be confined in the optical trap. Those in the upper hyperfine states will leave 
the trap by spin-flip scattering. Since spin-flip scattering is strong in 23 Na, only the high- 
field seeking stretched state |2, — 2) exhibits reasonable stability, experiments with more 
complex spinor condensate do not seem to be possible l36ll . On the other hand, optically 
trapped 87 Rb has proved to be a candidate for spin-2 Bose gas [ 37 ] with rich spin dynamics 
and magnetization conservation was also observed during the mixing [5]. In the case of 
85 Rb, the lowest multiplet has spin F = 2 and a negative s-wave scattering length in zero 
field. With the success to Bose condense 85 Rb in magnetic traps [7), it is conceivable that 
an F = 2 spinor condensate might be trapped optically in lower hyperfine states, provided 
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that the three particle losses when the field is reduced through the Feshbach resonance are 
not too large. 

Bose systems require that the total angular momentum of two colliding spin-2 particles 
is restricted to 0, 2, and 4. The effective low-energy Hamiltonian including the interaction 
energy describing binary collisions via the s-wave scattering can be generally expressed as 
(221 El [21 



H = /d 3 r|^(r)(-^ + C/(r))^(r)^ 



+-co^t(r)^(r)^(r)%(r) 
+id$J(r)*}(r) (F a ) ik (F a ) jt $ fc (r)tf,(r) 

+i5pa*J(r)*t(r) (2i;2j|00) <00|2fe;2Z> # fe (r)#,(r)J (13) 

where (00|2fc; 21) is the Clebsch-Gordan coefficient for combining two spin-2 particles with 
m-F = k and I into a spin singlet |0, 0). The parameters 

A-irh 2 4a 2 + 3a 4 

co = 



m 7 
Airh 2 04 — a 2 
m 7 



4vrn 2 3a 4 - 10a 2 + 7a 
5c 2 = (14) 

m 7 

describe the density-density interaction, spin-spin interaction, and formation of the singlet 
pair, respectively. The spinor £(r) with five components £(r) = (C+2, C+l> Co> C-1j C-2) 
normalized to unity, determines the average local spin as (F) = ^(r)F^(r), and F are 
the 5x5 spin-2 matrices which obey the same commutation relations [F a , F\,\ = ie a b c F c . 
In the mean-field approach the properties of a spinor condensate are determined by the 
spin-dependent energy functional 

K spin = (c 1 (F) 2 + c 2 \e\ 2 - P (F z )+q(Ftyn (15) 

where = 2( +2 C-2 — C+1C-1 + Co represents a singlet pair of identical spin-2 particles 
and is invariant under any rotation. The parameters p and q are related to the individual 
level shifts by 

p = ^(E +2 -E- 2 ) + ^(E- 1 -E +1 ) 

q = -LtE+z+E^ + ^iE^ + E^-^Eo (16) 
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The Breit-Rabi formula |34] in the case of 23 Na or 87 Rb (F = 2 is the upper hyperfine state 
with higher energy) gives 

E+2 = -l + \( 1 + x ) 
1 1 



E +1 = - 8+ - 2 V~l 

1 1 



+ X + X 



2 



EO = -77 + ^VT + 



x- 



1 1 



2 



E-i = -- + -\A - x + x 
o z 

E-2 = + 

In weak field, the quadratic Zeeman splitting is always negative, i.e., g = —jqX 2 + O (x 4 ) . 
In the case of 85 Rb (/ = 5/2) the lowest multiplet has spin F = 2. From the individual 
level shift 



11/4 

E+2 = \ l + -x + x 2 

+ - 12 2 V 3 



11/2 

= \l + -x + x 2 

+ 12 2 V 3 

E = — ■Jl + x 2 

u 12 2 v 



1 1 / 2 _ 

E_\ = \ 1 x + x 1 

12 2 V 3 



11/4 _ 
E - 2 = -12-2V 1 "3 X + X 

we easily see q is always positive for 85 Rb at small field, q = ^x 2 + O (x 4 ) . Unlike in the 
case of spin-1, the whole p-q plane is accessible experimentally for a spin-2 condensate. 

The ground state magnetization must be aligned with the external field, i.e. along z-axis, 
implying (F) 2 = (F z ) 2 in eq. ( fT3t . Minimization of the spin dependent energy functional 
using the similar Lagrange multiplier method leads to three possible phases, one more com- 
pared to the spin-1 case. These phases are characterized by a pair of parameters |(F)| and 
| G | describing the ferromagnetic order and the formation of singlet pairs, respectively. For 
convenience, we only consider the linear Zeeman shift p due to the magnetic field: 

• Polar/Anti-ferromagnetic phases 



2 V V 4ci - c 2 V 4ci — c 2 



2 V V 2( Cl -c 2 )' V 2( Cl -c 2 ) 

P : e w » (0,0, 1,0,0) 
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with energies C2 — p 2 / (4ci — C2) , C2 — p 2 /4 (c\ — c-i) and C2 respectively. Here #j 
are arbitrary phases for the corresponding components. These states are energetically 
degenerate in the absence of the external field with energy C2 and parameters (F) = 
and |0| = 1. 

• Ferromagnetic phases 

F : e i6 + 2 (1,0,0,0,0) 
F' : e ie+1 (0,1,0,0,0) 

with energies 4ci — 2p and c\ — p respectively. This phase has a non-vanishing 
parameter |(F)| = 1 indicating the ferromagnetic order and |0| =0. 

• Cyclic phase 




with energy — p 2 /Aci and cj> an arbitrary phase. This is a nonmagnetic phase which 
has no spin- 1 analog and was referred to as the cyclic state because of its close analog 
to the d-wave BCS superfluids. Both parameters are zero, (F) = and |0| =0. 

Recent experiments observed clear evidence of polar behaviour for F = 2 spinor con- 
densate of 87 Rb, and the slow dynamics of prepared cyclic ground states showed the F = 2 
state to be close to the cyclic phase ll26ll37l . However, the nature of the spinor conden- 
sate which depends on the s-wave scattering lengths for the total spins 0, 2, and 4, may be 
changed into other phases by an offset magnetic field. 

3 Outline of the homotopy theory of defects 

We sketch out the procedure which has been widely used in the study of topological defects 
in ordered media such as liquid crystals, superfiuid 3 He and heavy-fermion superconduc- 
tors. The explicit use of homotopy for topological classification of defects was made by 
some French ll38l l39l l40l l4~Tl and Russian authors ll42ll43ll . The results were well summa- 
rized in two review articles [39 44]. The central feature of the classification scheme of the 
defects emerges from examining the mappings of closed curves in physical space into the 
order-parameter space (OPS). 

The order parameter of a system has associated with it a group of transformations G. 
The set of all transformations in G that leave the reference order parameter / (chosen arbi- 
trarily but thereafter fixed) unchanged is known as the isotropy group H = {g £ G\gf = 
/}. The OPS can then be taken to be the space of cosets of H in G: M = G/H. In terms 
of broken symmetry, the fact that the ordering breaks the underlying symmetry is expressed 
in the fact that H is only a subgroup of G. The description that follows will be valid for 
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any group G that acts transitively on M, i.e., if fx and fi are possible values of the order 
parameter, then there is a transformation g in G which takes f\ into f^. = gfi- 

Homotopy groups of the order-parameter space describe physical defects [44]. The ra-th 
homotopy group ir n (M) of the space M consists of the equivalence classes of continuous 
maps from n-dimensional sphere S n to the space M. Two maps are equivalent if they are 
homotopic to one another. In three dimensional space, the first homotopy group, also called 
the fundamental group, 71"! (M) describes singular line defects and domain walls, which are 
non-singular defects. The second homotopy group ^(M) describes singular point defects 
and non-singular line defects. These can be calculated with the help of the fundamental 
theorem: Let G be a connected, simply connected continuous group and Hq be the set of 
points in H that can be connected to the identity by a continuous path lying entirely in H. 
Then we have the isomorphisms 



7Ti(Af) = H/Hq, 7T 2 (M) = TTiGffo). 



(17) 



For the theorem to hold, it is necessary that 7To(G) = 7Ti(G) = ^(G) = 0, meaning 
that G has only one connected piece, any loop in G can be shrunk continuously to a point, 
and G has a vanishing second homotopy group. While the second homotopy groups are 
always Abelian, the fundamental groups can either be Abelian (each element constitutes a 
conjugacy class), or non- Abelian (the line defects are characterized by the conjugacy classes 
instead of the elements). In Figure 1 we give a schematic description of the procedure for 
calculation of homotopy groups. 
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G ^H:g/=/^ H • ^^(G/H) ^^(G/H) 



Reference 
Spinor / 



OPS:G/H 



Identity / 

Figure 1: A schematic description of the procedure for calculation of homotopy 

groups. 

A ready example for illustrating the above procedure is the biaxial nematics, whose 
symmetry is that of a rectangular box (proper point group D2). If G is taken to be SO(3) 
then the isotropy subgroup H is the four-element group consisting of the identity and 180° 
rotations about three mutually perpendicular axes {D2). Order parameter space is thus 
identified as M = SO(3)/L>2- If. however, we take G to be SU(2), the universal covering 
group of SO (3), then H is expanded to the non- Abelian quaternion group Q (known as the 
lift or double group) with eight elements 



Q = {±1, ±io- x ,±ia y , ±icr z }. 



(18) 



The natural representation for the order parameter space of a biaxial nematic turns out to be 
M = SU(2) IQ. Since it is a discrete subgroup of SU(2), H/H = H. Thus m(M) = Q, 
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and 7T2(M) = 0. There are no stable point defects in biaxial nematics and the line defects 
are characterized by five conjugacy classes of group Q 



Co = {1},C = {-1}, 

C x = {±ia x }, C y = {±ia y }, C z = {±ia z }. 



(19) 



The class Co contains removable trivial defects; Co contains defects in which the object 
rotates about 360° as the defect line is encircled; the other three classes contain defects in 
which the rotation is through 180° about each of the three distinct symmetry axes. The 
defects here are non-commutative, providing an example with non-Abelian fundamental 
group. 

Another illustrative example is the dipole-free ^4-phase of 3 He, which affords an un- 
usual example of a case where C must be bigger than SO(3). The order parameter is the 
product of an arbitrary unit 3-vector n and a complex 3-vector of the form u + iv, where 
u and v are an orthonormal pair. The orientations of n and u + iv are uncoupled. Take 
the reference order parameter to be Ay = zi(xj + iyj), the group C can be taken to be 
the direct product of SO(3) with itself: G = SO(3) x SO(3), elements of G consisting 
of pairs (R, R') of distinct rotations. The isotropy group H consists of elements of the 
form (R(z,9), 1) and (R(u, ir), R(z, ir)) for any axis u in the x-y plane. To construct a 
simply connected G, we must replace each SO (3) by SU(2). Determining the lift of H 
from SO(3) x SO(3) to its covering group SU(2) x SU(2), we find the isotropy group 
consists of 4 pieces {Hq, gHo : g 2 Ho, g 3 Ho} with the connected component of the identity 
Hq = {(u(z,9), 1)} and g = (u(x,ir) ,u(z,ir)). In this article, the notations R and u 
represent the rotations in SO(3) and SU(2), respectively. The fundamental group is thus 
isomorphic to the cyclic group of order 4, tt\{M) = and ^(M) = Z. 

For the spinor condensate it seems natural to identify the underlying symmetry group as 
U(l) x SO(3), the groups in the direct product representing the gauge and spin degrees of 
freedom respectively. This group is not simply connected, i.e., 7Ti(U(l) x SO(3)) / 0. To 
apply the theorem, however, it is again essential that one chooses the group G to be simply 
connected. We proceed by specifying the symmetry group as its universal covering group 
R x SU(2), with the group of real numbers R representing any translation 6 £ (— oo, +oo) 
in the phase of the condensate. For F = 1, we use the 3D representation of the group SU(2) 
in order to obtain the isotropy group, e.g., a rotation u(z, a) around axis z by angle a takes 
the form of a diagonal matrix Diag {e~ ia , 1, e ia ) , a rotation u(y, (3) around axis y by angle 
[3 takes the form of 



even representation of SU(2), though we know (and should always bear in mind) that 

u(z, a + 2ir) = —u(z, a) while u(z, a + 4-7r) = u(z, a). 




, l,e ia ) in this 
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4 Calculation of the homotopy groups 



There are two possible ground states in F = 1 case. For the ferromagnetic state, the 
isotropy group H is constructed by the set of transformations which leave the reference 
order parameter (1, 0, 0) T invariant. From the degenerate family of the ground state spinor 
eq.© we know immediately that the angles should satisfy 

(3 = 0, 9 - a - 7 = 2mr (20) 

with n an integer. The elements in group H are the combination of a translational part and 
a rotational part H = {(6>, u(z, 9)) , (9, u(z, 9 + 2ir))} = {(9, ±u(z, 9))}. Evidently this 
group includes two disconnected components — the connected component of the identity 
Hq = {(9,u(z,9))} is isomorphic to R. The group H/Hq is isomorphic to the integers 
modulo 2, i.e., Z 2 . The second homotopy group 7r 2 is trivial and we arrive at the same result 
as that in Ref. 021 

vri(M) = Z 2 , vr 2 (M) = 0, (spin-1 FM state). (21) 

A ferromagnetic spin- 1 condensate may have therefore only singular vortices with winding 
number one while the point-like defects are topologically unstable. Alternatively we may 
take the symmetry group as SU(2) because we can produce all possible gauge transforma- 
tions by absorbing 9 into the Euler angle 7. The isotropy group is discrete and isomorphic 
to Z 2 , which gives exactly the same result. 

The polar state emerges if the atoms in the condensate interact anti-ferromagnetically. 
In the ground state eq. ©, the reference parameter (0, 1, 0) T is left invariant for just those 
elements with 

p = 0,9 = 2mr or = vr, 9 = (2n + 1) vr. (22) 

Thus the isotropy group H includes now the transformations in which both the rotation 
and the translation leave the spinor unchanged, and those in which the rotation takes the 
reference spinor (0, 1,0) T to (0, — 1, 0) T and the translation takes it back, i.e., a tt rota- 
tion about arbitrary axis perpendicular to z combined with a tt translation in (or any 
odd multiples of tt). The latter invariance is identical to the Ising gauge symmetry em- 
phasized in eq. (14) of Ref. [20|. The full isotropy group is the union of these two 
sets, H = {(2n7r, u(z, a)) , ((2n + 1)tt, gu(z, a))} where g = u(u,tt). There are in- 
finitely many discrete components in H, while the connected component of the identity 
Hq = {(0, u(z, a))} is isomorphic to U(l). The elements with an even translational parity 
are of the form (2mr, I)Hq, and those with an odd parity are of the form ((2n + 1)tt, g) Hq. 
The group H/Hq is therefore isomorphic to the group of integers Z through the isomor- 
phism ((2n + g 3 \ Hq i— > 2n + j for j = 0, 1. We recover the conclusion that line and 
point defects in spin- 1 polar state can be classified by integer winding numbers, 

7Ti(M) = Z, tt 2 (M) = Z, (spin-1 Polar state). (23) 

Thus the Z2 term does not appear in the homotopy group. We argue that the identification 
of the OPS in Ref. [ 19 20 1 is also incorrect (see below). Physically there are indeed infinite 
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number of line defects corresponding to integer and half-integer vortices (eq. (27) in Ref. 
[20 1). On the other hand, it is the Ising symmetry that leads to half- vortices (j = 1), which 
have been shown to be the unique linear defects in polar condensate in addition to the usual 
integer vortices (j = 0) [18]. If we move around a closed path in the condensate we note 
that when we return to the starting point the angle 9 has changed by some amount. If we 
define the change in this angle divided by 2ir to be the winding number, we see from the 
elements of H/H that the winding number can be either an integer n or a half-integer 
n + 1/2. 



5 Spin-2 Bose condensate 

We next apply the same approach to the BEC of spin-2 bosons. The defects which may be 
created in spin-2 condensate exhibit even more elaborate structures due to quantum corre- 
lations among bosons. For F = 2 we have to use the 5D representation of SU{2), e.g., the 
rotation u(z, a) is represented by matrix Diag (e 
the form of 



-2ia 



) and u(y,/3) takes 



4 
COS ^ 

sin (5 cos 2 | 
^ sin 2 
sin P sin 2 ^ 

■ 4 



— sin P cos 2 | 

c os /3+cos 20 
o 



2 

^ sin 2/3 

cos — cos 20 
2 

sin P sin 2 ^ 



^ sin 2 P 
-^sin2/? 

1+3 cos 20 
4 

^ sin 2/3 



4 

' 2 /3 



^sin< 



- sin P sin 2 | 

cos 0— cos 20 
2 

sin 2/3 

s /3+cos 20 
2 

./?COS 2 f 



sm 



4 



cos / 

sin 



sin P sin 2 ^ 

P 

sin p cos 2 ^ 
cos 4 f 



^sin 2 



/ 



The calculations of the degenerate family of the ground state spinors and the corresponding 
homotopy groups are straightforward and some results have been reported in [47 1. Here we 
pick up some interesting features in our results, focusing on the symmetry properties of the 
defects in comparison with those in other ordered media. We first consider the defects in 
the absence of an external field and the effect of magnetic field will be discussed later. 

We start with the case of the ferromagnetic state F. Equating the general expression for 
the ground state spinor 



Ad-2-y) 



( 



\ 



e 

—ia 



2ia 



cos 4 f \ 



sin p cos 2 f 
^ sin 2 P 
e ia sin P sin 2 j 
e 2ia sin 4 f 



(24) 



/ 



with the reference spinor (1, 0, 0, 0, 0) leads to the requirement for the isotropy group H 



P = 0, 9 - 2a - 27 = 2nir. 



(25) 



We see that taking n = 0, 1, 2, 3 is enough for all possible transformations, with the trans- 
lational part being arbitrary and the rotational part containing the rotations around axis z 
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by 6/2,9/2 + it, 9/2 + 2it,9/2 + 3ir respectively. Hence the group H is composed of 
four pieces H = {(9, u(z, nir + 9/2))}. Here it is important to show that the four compo- 
nents are not connected: there does not exist a continuous path in H which connects one 
component to another, though the rotational parts themselves are connected. The connected 
component of the identity H Q = {(9, u(z,9/2))} is again isomorphic to R. If we define an 
element g of the group R x SU(2) by (0, u(z, it)), we see that the quotient group H/Hq has 
the same structure as the cyclic group of order 4, i.e., {e, g, g 2 ,g 3 } and we conclude that 



tti(M) = Zi, 7T 2 (M) = 0, (spin-2 F state). 



(26) 



It is interesting to check how the group Z4 characterizes vortices for state F. In spin-1 
case there is only one topologically stable line defect, that is, a vortex with winding number 
one. Equation d26t shows that there are three stable vortices for spin-2 condensates. We 
can set 9 — 27 = 2rrup, —a = imp, = nt in the ground state for F state, Eq. d24l . 
which leads to a family of spinor states parametrized by a parameter t between and 1. 
Here m > is an integer, tp is the azimuthal angle. When t evolves from to 1, the 
Amp vortex state ((t = 0) = (e* 4mip , 0, 0, 0, 0) T evolves continuously to the vortex free 
state C(t = 1) = (0,0,0,0, 1) T . This shows that vortices with winding number 4m are 
topologically unstable. Similarly, by multiplying factors e lkip (k = 1,2,3) one obtains the 
following correspondences 



(27) 



i.e., the vortices with winding numbers 4m + k may evolve into vortices with winding 
numbers k, respectively. There are thus three classes of topologically stable line defects. 
Together with the uniform state, they form the fundamental group Z4. Non-trivial vortices 
are those in which the reference spinor rotates through 180°, 360° or 540° about the z-axis 
when the defect line is circulated. Straightforwardly for ferromagnetic condensates with 
spin F, the fundamental group tti(M) = Z2F characterizes (2F — 1) classes of stable line 
defects. 

Spin variations in the ferromagnetic states in general lead to superflows |3j|25|. To 
illustrate the coreless (or Skyrmion) vortices in spin-2 case, we set 9 — 2r = 2ip, a = in 
the spinor degenerate family d24l and consider the condensate 



C(r) 



4 

COS 

sin cos 2 f 
e i2 ^sin 2 



\ 



P 

e i3< ? sin sin 2 f 



V 



sin' 6 



,4 £ 
2 



where = 0(r) is an increasing function of r starting from = 
velocity does not depend on z and it is cylindrically symmetric 



(28) 



at r = O.The superfluid 



h . 



2 cos 0V(p] 



2h . 



(29) 
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i.e., the coreless vortex may exist in the spin-2 case, with only the velocity doubling its 
value compared to the spin-1 case 0. The velocity vanishes instead of diverging at r = 
because (3(0) = 0. This is called a coreless vortex. For a Mermin-Ho vortex [45), the 
bending angle (3 must be tt /2 at the boundary of the condensate, while for an Anderson- 
Toulouse [46 1 vortex (3 must be tt, i.e. 

(3(R) = vr/2, for Mermin-Ho 

(3(R) = tt, for Anderson-Toulouse. (30) 



6 Non-Abelian homotopy groups 

Media with non-Abelian fundamental groups are especially interesting from the topological 
point of view. The only illustrative example in ordered media so far have been biaxial 
nematic liquid crystals [48 1. Their multiplication table has been verified experimentally 

ED. 

We have found that the cyclic state C provides another physically realistic example 
in which the fundamental group is non-commutative. A rotation and a gauge transforma- 

rp 

tion of the reference spinor | in zero field produce the following 

degenerate family 



■ -2iat / cog 4 §_ e i(p-2iT _|_ \/§ , 



4 P -i(/>+2iT 



< -'" I c-os ■ -r" ' -"■ + ^ sin (3 — sin 
e~ ia sin (3 (cos 2 f e^" 2iT - V3 cos (3 + sin 2 i e -^+ 2iT 



^ sin 2 (3e^- 2iT + ^(1 + 3 cos 2(3) - & sin 2 (3e 
e ia sin (3 (sin 2 ^~ 2iT + V3cos/3 + cos 2 f e "^ +2 " 

r + ^ sin 2^_ cos 4| e -^+2i 



) \ 



V 



-i<f>+2ir 



e 2ia ( sin 4 § e^" 2 " 



/ 



The reference spinor is left invariant by the elements of three sets characterized by the 
translations in the phase of the condensate 9 : 



For 9 = 2n7r, one must have (3 = 0, a + r = m7r, or (3 = tt, a — r = — + § + run; 
For # = + 2mr, one must have (3 = + t = — ^ + mix, a — r = —<f> + mV; 
For 9 = ^ + 2n7r, one must have /3 = : |,a + r = ^ + rmr, a — t = —cf> + rn'ix. 



Here m and m! are integers satisfying m + m' = odd. For all possible transformations 
we need take the values m, m' = 0, 1, 2, 3 so that there are eight possibilities 



m 


= 0,m 


= 1,3 


m 


= 1, mf 


= 0,2 


m 


= 2,m 


= 1,3 


m 


= 3, m' 


= 0,2 
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This gives the isotropy group 

H = {±I,±a,±b,±c, (31) 
±d,±e,±f, ±g, 
±d 2 ,±e 2 ,±f 2 ,±g 2 }. 

The spin rotations a = u(z, tt), b = u(y, n)u(z, <fi + vr/2) and c = ba satisfy a 2 = b 2 = 
c 2 = —I, while d = u(z, 7r/4 + (p/2)u(y, n/2)u(z, 7r/4 — (f>/2), e = —da, f = —ad and 
g = —ada satisfy d 3 = e 3 = f 3 = g 3 = —I. Each element in the first, second, third row is 
associated with an additional phase change 2mr, 2tt/3 + 2nir, Air/3 + 2wk respectively. It 
is a discrete group, and Hq consists of the identity (0, I) only. The fundamental theorems 
identify that 

tti(M) = H, vr 2 (M) = 0, (spin-2 C state). (32) 

The elements in the fundamental group are non-commutative, for example ab = —c / ba. 

The criterion for the topological equivalence of defects applies in the most general case 
in terms of conjugacy classes of the fundamental group. Two line defects are topologically 
equivalent if and only if they are characterized by the same conjugacy class. Defects can 
still be labelled by the elements of the first homotopy group, but if these elements belong 
to the same conjugacy class, corresponding defects can be continuously transformed to one 
another. However, if they belong to different conjugacy classes this is not possible. It is 
thus necessary to classify the group into the following conjugacy classes: 

C (n) = {/}„, Cb(n) = {-/}„, C 2 (n) = {±a,±b,±c} n , 
C 3 (n + l/3) = {d,e,f,g} n+1 / 3 , C 3 (n + 1/3) = {-d, -e, -/, -g} n +y 3 , 
C|(n + 2/3) = {d 2 ,e 2 J 2 ,g 2 } n+2/3 , C 2 (n + 2/3) = {-d 2 , -e 2 , -f 2 , -g 2 } n+2/3 

with the subscripts standing for the winding numbers of the defects. Physically this indi- 
cates the feasibility of creating not only vortices with any integer winding number but also 
fractional quantum vortices. The class Cb(n) describes defects in which the phase of the 
spinor is changed by 2im as the defect line is encircled. Note that only Cb(0) corresponds 
to trivial defects. In the case of Co(n) phase change of 27rn is accompanied by a 360° rota- 
tion about z-axis. The element a with winding number n in the class C2(n) depicts a defect 
in which the spinor rotates through 180° about the z-axis and changes phase by 2im as the 
line is encircled. The multiplication table of conjugacy classes is shown in table 1. Only 
half of table is shown because the class multiplication is commutative. Winding numbers 
have been omitted for clarity. When two classes are multiplied the winding number of the 
resulting class is the sum of the individual winding numbers. It shows that, for example, 
when we combine defect C 2 {n) with C 2 {— n) they can either annihilate each other (Cq(0)) 
or form defect Cb(0) or C 2 (0), the result depending on how they are brought together. In- 
teresting features of this non-Abelian fundamental group include the topological instability 
of the defects and their interaction, i.e., entanglement when two of them are brought to cross 
with each other [44]. 
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Table 1: The multiplication table of the conjugacy classes of the cyclic phase. 





Co 


c 2 


c 3 


c 3 


c 3 2 ci 


Co 


C 










c 2 


c 2 


QC + 6Co + 4C 2 








c 3 


c 3 


3{C 3 + C~ 3 ) 








c 3 


c 3 


3 (C 3 + C 3 ) 


Ci + 3C| 


3C| + C| 




ci 


ci 


3 (Ci + cf\ 


4C^ + 2C 2 


4C + 2C 2 


3C 3 + C~ 3 


ci 


ci 


^ (ci + cl) 


4C + 2C 2 


4Cb + 2C 2 


C 3 + 3C 3 3C 3 + C 3 



The defects can be further grouped into classes, which form an Abelian group isomor- 
phic to the first homology group of the order parameter space 1 50 1. This coarser classifica- 
tion is more general than the homotopic one because two defects are considered equivalent 
also, if they can be transformed into each other via a catalyzation process consisting of 
splitting a line singularity into two and recombining them beyond a third one. All elements 
labelled by elements of the commutator subgroup D of tv\ (M) can be catalyzed away by 
this procedure. D is generated by the commutators 5t5~ 1 t~ 1 of all pairs of elements 
5, t 6 tt\(M). The elements of tti(M)/D are unions of conjugacy classes. In our case D 
is the union of the conjugacy classes with winding number zero, D = Co(0)UCo(0)UC 2 (0) 
and the first homology group is 

TTi (M)/D = [C U C~ U C 2 , C 3 U C~ 3 , Cl U Cl] (33) 

The homology theory assembles the conjugacy classes further into three sets for each n, in 
which the defects are labeled by the winding numbers n,n + l/ 3,n + 2/3 respectively K71 . 
Two defects in the same conjugacy classes can be continuously converted into one another 
by local surgery, while two defects in the same homology class can be deformed into one 
another by the catalyzation procedure. 

7 Order Parameter Spaces 

Like the quaternion group Q for biaxial nematics, the fundamental group d3*Tl is the lift 
of a point group in R x SU(2). To find the remaining discrete symmetry group for the 
cyclic state, and, in addition, to clarify the controversial identification of the OPS for spin-1 
case, in the remaining of this paper we turn to describe the system in terms of rotations in 
SO(3), e.g., two elements ±u(z, a) in SU(2) are mapped into one R(z, a) in SO(3) with 
R(z, a + 2ir) = R(z, a). 

The OPS for F = 1 polar state was identified as U(l) x S 2 in Refs. EEl. An extra Z 2 
symmetry was claimed in Ref. [ 19 ] so the author concluded the OPS as U(l) x S 2 /Z 2 . Here 
we show that previous studies are incorrect. Taking the group G as U(1) G x SO(3) s where 
the subscripts stand for the gauge and spin symmetries respectively, we see that the isotropy 
group H consists of two separate parts, {(e*°, R(z, a)) \ and { (e t7r , R(y, ir)R(z, a)) }. 
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The rotations in the first part constitute the group SO (2), while the elements in the second 
part are just those in the group 0(2) but not in S0(2) with determinants —1. The combi- 
nation of these two parts gives the full isotropy group as 0(2) where both gauge and spin 
symmetries are involved. The OPS is the quotient G/H = (U(1) G x SO(3) 5 ) /0(2) G+S 
and here it is not possible to apply the fundamental theorem for G is not any more simply 
connected. One may wonder if we can factorize the OPS further as 

G/H = (U(l) G xSO(3) s )/0(2) G+s 

= (U(l) G xSO(3) s )/(SO(2)xZ 2 ) G+5 
= U(l) x S 2 /Z 2 

However it is incorrect because though in 3 dimensional space we have 0(3) = SO (3) x Z 2 
but it is not true in 2 dimensional case, i.e. 0(2) / S0(2) x Z 2 . The spin and gauge 
symmetries are broken in a connected fashion just as in the system of 3 He 03 ED- Table 2 
summarizes our result in comparison with the previous studies. 



Table 2: Comparison of the OPS and fundamental groups for spin-1 polar condensate 
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vri(M) 


Ho, Stoof, etc. 


U(l) 


x S 2 
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Zhou 


U(l) 


x S 2 /Z 2 


Z x Z 2 


This paper 


(U(l 


) xSO(3))/0(2) 
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Figure 2: Symmetries of the defects in biaxial nematics (D 2 ) and cyclic state C in 
spin-2 condensate (T). The dot at the center of the rectangle stands for axis z. The 
dashed lines represent 2-fold axes, except that with a triangle for 3-fold axis. 

For the ferromagnetic state the group H may be obtained if one notices that the 2ir 
difference in the rotational angle does not give another component as it did in the case of 
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SU(2). We have H = { (e i0 , R(z, 9)) } which is isomorphic to U(1) G+S . This means that 
there is a remaining symmetry U(l) in the symmetry broken system. The OPS is thus 
factorized as (U(1) G x SO(3) s ) /U(l) G+5 = SO(3) s+G . 

The discrete symmetry group of defects in the spin-2 cyclic state C can be shown to be 
isomorphic to the tetrahedral group T. We continue to represent G as U(1) G x SO(3) s . 
The isotropy group d3TI > is shrunk to a group of 12 elements if one understands the rotation 
in the sense of SO(3) (i.e., a = R(z, it)), 

H = {I, a, b, c, ed, ee, ef, eg, e 2 d 2 ,e 2 e 2 ,e 2 f 2 , e 2 g 2 }, (34) 

where e = exp(2-7rz/3) comes from the gauge transformation and ed, for instance, is an 
abbreviation for the element {e,d). Three 2-fold rotational axes are z and 2 lines in xy 
plane perpendicular to each other (which lie on axes x and y if we choose the arbitrary phase 
4> = 7r/2). The elements ed, ee, ef, eg are four 3-fold axes. The symmetries remaining in 
the symmetry broken states for biaxial nematics and spin-2 cyclic state are shown in Figure 
2. The OPS for state C can be identified as (U(1) G x SO(3) s ) /T G+S . 

It should be noted that an applied magnetic field B changes the defect structure severely 
by reducing the degenerate family of the spinor. We take again the cyclic state as an exam- 
ple. The symmetry group in this case is U(l) x SO(2) because the magnetic field chooses 
its direction automatically as the quantization axis. From the spinor 





y/2 - 2p2 e ^l+<«/2 





where <pi : 2 are two arbitrary phases and p ~ B, we easily see the possibility to create 
vortices in any of the three nonzero components with winding number for i— th component 
rii confined by m + = 2ri3. 



8 Summary 

Our main findings are summarized in Table 3. We have determined the nature of the topo- 
logical defects in spin-1 and spin-2 condensates. The order parameter spaces are identified 
as the spaces of the coset of the isotropy group H in the transformation group G. Topo- 
logical^ stable vortices with winding numbers larger than unity may be created in the 
ferromagnetic state for condensates with F > 1, up to the value (2F — 1). The line defects 
in the spin-2 cyclic state C exhibit non-commutative features, resulting e.g. in line de- 
fects with winding numbers of 1/3 and its multiples. It also turns out that in the zero field 
U(l) x SO(3) does not act transitively on the order-parameter space of the polar phase and 
thus the defect structure remains unsolved. 
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Table 3: Main results on calculation of the OPS and homotopy groups 
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Spin-2 F 
Spin-2 F' 
Spin-2 C 



Spin-1 FM 
Spin-1 AFM 
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Z 2 

z z 

Z 4 

z 2 

H eq.(EB 
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